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Abstract. The hydrodynamic model is used to describe the single-particle 
spectra and elliptic flow of hadrons at RHIC and to predict the emission angle 
dependence of HBT correlations at RHIC and LHC energies. 



1. Introduction and Summary 

One of the key signatures of the formation of thcrmalizcd hot and dense matter, 
such as a quark-gluon plasma (QGP), in relativistic heavy-ion collisions is collective 
transverse flow, driven by the transverse pressure gradients between the center of the 
hot collision fireball and the surrounding vacuum. The space-time evolution of this 
collective flow can be described by relativistic hydrodynamics. In the limit of zero 
mean free path (i.e. when microscopic scattering rates exceed any macroscopically 
relevant time scale) the system can be described as an ideal fluid neglecting transport 
effects such as viscosity and heat conduction. The flow observed in the final state via 
analysis of the single-particle spectra of emitted hadrons can then be directly related 
to the initial conditions of the collision fireball when it first thermalized, and to the 
equation of state p(e, ns) of the expanding matter. 

In non-central collisions the initial conditions are anisotropic, due to the spatial 
deformation of the nuclear overlap region in the transverse plane, resulting in 
anisotropic pressure gradients which ultimately lead to anisotropies in the final 
momentum distributions. The momentum anisotropies can be quantified by the 
coefficients v n of an azimuthal Fourier decomposition of the transverse momentum 
distributions. These coefficients, in particular the elliptic flow coefficient V2, grow 
only as long as anisotropies in the spatial pressure gradients persist. Since the 
reaction zone quickly evolves from its initial spatially deformed state towards a more 
cylindrical shape without significant anisotropies, the growth of vi is restricted to 
the earliest stages of the fireball expansion £Q. Even in the absence of interactions 
among the fireball constituents the spatial anisotropy which could drive the build-up 
of elliptic flow decays quickly, due to radial free-streaming of the matter 0. Delayed 
thermalization therefore manifests itself in reduced elliptic flow since vi then has 
to build on smaller anisotropies in the pressure gradients. If thermalization happens 
quickly, elliptic flow builds up earlier and reaches larger values. The resulting stronger 
flow into the reaction plane [5J destroys the spatial deformation even faster than free- 
streaming, resulting in an early saturation of the momentum anisotropy after about 5 
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fm/c in semiperipheral Au+Au collisions at RHIC For a given initial deformation 
(i.e. a fixed impact parameter) the saturation value of v 2 has been shown to be a 
monotonic function of the microscopic scattering cross section [1], reaching a strict 
upper limit given by ideal fluid hydrodynamics in the limit of zero mean free path. 
The fact that the elliptic flow measured at RHIC almost exhausts this hydrodynamic 
limit (for an extensive list of references I refer to the recent review 0) is at this 
point the strongest available evidence for early thermalization (at Ttherm < lfm/c [Hj) 
and, together with an estimated energy density > 10GeV/fm 3 at this early time 0, 
strongly suggests the formation of a well-thermalized QGP with significant (5-7 fm/c) 
lifetime 0. 

We have demonstrated elsewhere |H] that the bulk of the observed hadrons 
from Au+Au collisions at RHIC (i.e. all mesons in the region pr < 1.5 — 2GeV/c 
and all baryons at pt $ 2 — 2.5GeV/c) behave hydronymically, i.e. that their single 
particle spectra and momentum anisotropies in central and semicentral collisions with 
impact parameter b < 7 fm can all be quantitatively reproduced by relativistic ideal 
fluid hydrodynamics with a common kinetic decoupling temperature. In Section 
2 I supplement this finding by pointing out that the same holds true even for the 
heavy f2 and Q hyperons, in contrast to earlier expectations based on their presumed 
weaker coupling to the expanding pion fluid during the late hadronic expansion stage 
[7]. It has also been noted El that the hydrodynamic model does much more 
poorly in its description of the two-particle Bose-Einstein (Hanbury Brown Twiss 
or HBT) correlations, by overpredicting the longitudinal and outward HBT radii 
while underpredicting the sideward radius and its transverse momentum dependence 
(the "RHIC HBT puzzle")- In Section 3 I show how HBT studies as function of 
the azimuthal emission angle can constrain the space-time evolution of the spatial 
deformation in non-central collisions, thereby supplementing the information from the 
measured momentum-space anisotropies and constraining the total duration of the 
expansion phase. Preliminary RHIC data from 2-pion correlations E3j agree with 
the hydrodynamic prediction jll| that at kinetic freeze-out the fireball is (on average) 
still slightly out-of-plane elongated, although much less so than initially. This is 
expected to change at LHC energies where the much higher initial energy density 
allows the elliptic flow to act for a significantly longer time period until decoupling, 
leading to pion emission from an (on average) in-plane elongated source. The HBT 
signatures of such an in-plane elongation are discussed in Section 3. 

2. Radial flow of 51 hyperons 

When shortly after the beginning of the experimental program at RHIC in 2000 the 
first data on the total charged particle multiplicity in central Au+Au became available 
|12| . Peter Kolb, Pasi Huovinen and I used this to constrain the initial conditions 
for the hydrodynamic model in central conditions and to predict hadron spectra 
and their anisotropies for non-central collisions Some of these predictions are 

shown in Fig. The three sets of curves correspond, from top to bottom, to central 
(6<5.4fm), semiperipheral (5.4 fm < b < 9.9 fm) and peripheral (9.9 fm < b < 13.5 fm) 
Au+Au collisions at y/s — 130 A GeV. At that time no spectra were available yet to 
fix the freeze-out temperature, so we estimated it to probably lie between 120 and 140 
MeV (solid and dashed curves, respectively). We later learned [5] that Tf = 130 MeV 
fits the pion and proton data from central collisions best. For the hyperons 
in the lower right panel we also included a prediction for freeze-out directly after 
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hadronization at T c = 165 MeV, motivated by the prevailing prejudice \7\ that, due to 
their lack of resonant scattering with pions, they would not be able to efficiently pick 
up any additional transverse flow generated in the hadronic phase. This expectation 




Figure 1. Hydrodynamically predicted hadron spectra from Ref. 1131 . together 
with preliminary f2+Q spectra from the STAR collaboration [141 . See text for 
details. 



was supported by observations by the WA97 Collaboration at the SPS JS] which 
indicated a smaller slope parameter for and than for other hadrons. It was only 
recently realized that this evidence was at least partially distorted by fitting a curved 
Q spectrum in a limited wit region by a simple exponential. Recent NA49 data for £1 
and Cl at smaller hit show evidence that the spectra bend over at low mx and seem to 
be consistent with much larger radial flow than previously thought, almost as strong 
as that experienced by pions and nucleons |16| . 

Figure ^ shows that the same is true at RHIC: The fi spectra from Au+Au 
collisions at t/s = 130 and 200AGeV [2J are consistent with the hydrodynamic 
predictions if they are assumed to freeze out at the same decoupling temperature 
T{ ~ 130 MeV as pions, kaons and (anti)protons, while they are much too flat to 
be compatible with freeze-out already at T c , i.e. directly after hadronization. It is 
impossible to generate enough radial flow for the f2s already before hadronization, 
without destroying the agreement of the model with all other hadron spectra and/or 
lowering the hadronization temperature significantly below 165 MeV (which would 
contradict lattice QCD results). We conclude that the hyperons seem to be much 
more strongly coupled to the expanding pion fluid during the late hadronic stage 
than previously thought. It will be interesting to explore the microscopic mechanisms 
responsible for such a strong coupling. 
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3. Emission angle dependent HBT interferometry 

The observational constraints on dynamical models for the evolution of the heavy-ion 
reaction zone can be further tightened by a simultaneous analysis of the single-particle 
momentum spectra with 2-particle momentum correlations. The latter provide access 
to the space-time structure of the source at kinetic freeze-out. More specifically, in 
the case of Bose-Einstein correlations between identical bosons, the Gaussian widths 
Rlj (K) of the correlation function C(q, K) = 1 + \(K) exp [- E 4 , i=0 , s ,i HQjRi,j ( K )] 
are linear combinations of the components of the spatial correlation tensor S^ V {K) = 
[(Xfj,x u ) — (x^) {x u )\ (K) which characterize the space-time widths of the source 
emission function S(x 7 K) at freeze-out. For the general case these relations can be 
found in |17l I18| : for longitudinally boost-invariant sources and pions at midrapidity 
Y = they simplify to 

R 2 (3>) =±(S xx + S yy )-l(S xx -S yy )cos2$-S xy sin2<f> (1) 

R 2 ($) = \{S XX + S yy ) + \{S XX - S m )zos2§ + S xy sm2$ 

-2i3r(S tx cos$ + St y smty+faStt (2) 
R 2 os ($) = S xy cos 2$ - ±(S XX - S* TO )sin2$ 

+ p T {Stx sin 2$ -S ty cos$) (3) 

R?m =s zz . (4) 

For simplicity the dependence of the spatial correlation tensor S^ on the pair 
momentum K (in particular its implicit dependence on the azimuthal emission angle 
$ of if relative to the reaction plane) has been dropped, and only the explicit 
^-dependence arising from the azimuthal rotation between the outward direction 
x ou t || Kt and the reaction plane coordinate x \\ b is exhibited. For small transverse 
pair momenta Kt this explicit ^-dependence dominates ^Hj- The source widths S^ 
are specified in coordinates (x, y) parallel and perpendicular to the reaction plane 
while the correlation function is measured in (cc ou t, a^side) coordinates defined by the 
directions parallel and perpendicular to Kt- Pt is the transverse pair velocity. The 
terms R? ol and R 2 sl vanish due to the assumed longitudinal boost-invariance. 

The orientation of the reaction plane in the laboratory system can be determined 
event-by-event from the elliptic flow. The correlation function can thus be measured 
as a function of the azimuthal emission angle $ relative to the reaction plane. For 
non-central collisions the spatial deformation of the source in the transverse plane 
leads to oscillations as a function of $. In practice the finite angular bin size and 
reaction plane resolution smear these azimuthal oscillations and reduce their measured 
amplitudes, but this can be fully corrected for in a model-independent way |18| . 
General symmetries ensure that the HBT radii extracted from a Gaussian fit to the 
thus corrected correlation function exhibit the following general ^-dependences |18j : 

R 2 M =<o + 2E„= 2 ,4 > 6,...^,nCos(n$), 

Rim =<o + 2E„ = 2,4,6,... i? o,™COs(n<f), 

Rim = 2 E^ja C sin(n$), 

Rf^) - Rlo + 2En=2A^..RtnCO S (n®). (5) 

For sources without longitudinal boost-invariance R a i and R 2 ol oscillate around zero 
with cosines and sines, respectively, of odd multiples of $ |18l IB?] . The relations 
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continue to hold true for correlation measurements made in a finite symmetric window 
around Y = 0. 
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Figure 2. Contour plots of and horizontal (solid) and vertical (dashed) cuts 
through the .fC-integrated pion emission function in the transverse plane for RHIC- 
and LHC-typical initial conditions, for Au+Au collisions at 6 = 7fm 1111 . 



Figure [21 shows the momentum-integrated hydrodynamic emission function for 
RHIC- and LHC-type initial conditions. The curves labelled "RHIC 1" correspond to 
the same initial conditions which led to a successful description of the single-particle 
momentum spectra and elliptic flow in 130AGeV Au+Au collisions. One sees that, 
when integrated over all momenta and all emission times, this source is still out-of- 
plane elongated, although only very slightly. "IPES" labels an in-plane-elongated 
source obtained by starting the hydrodynamics earlier at a higher initial temperature 
of 2GeV In this case the source lives longer, grows much bigger and changes its 
initial out-of-plane elongation into a significant in-plane elongation at freeze-out. 

The vertical and horizontal cuts in the right panel of Fig.|2show that both sources 
are very "opaque" in that most particles are emitted from the (deformed) surface. 
A differential analysis shows that only pions with small transverse momenta 
come from the interior of the source; they are thus the only ones which probe the 
overall spatial deformation of the source in the transverse plane. Pions with non-zero 
transverse momenta, on the other hand, come from thin slivers near the edge of the 
fireball in the emission direction (for details see JJ); they do not directly probe the 
overall shape of the source. 

Figure [3] shows the azimuthal oscillations of the four non- vanishing HBT radius 
parameters for several values of Kt, calculated from these emission functions. The 
opposite spatial deformation of the RHIC 1 and IPES sources is reflected in the 
opposite signs of the oscillations for R^, R% and R^ s at Kt = 0- For example, at 
RHIC energies R 2 S oscillates downward, implying a larger sideward radius when viewed 
from the x direction (i.e. within the reaction plane) than from the y-direction (i.e. 
perpendicular to the reaction plane). For the in-plane elongated source (IPES) some 
of the oscillation amplitudes change sign at larger transverse momenta. This change of 
sign originates from an intricate interplay between geometric, dynamical and temporal 
aspects of the source at freeze-out \H\ and reinforces the statement that at ^ 
the HBT correlations do not probe the entire source and its spatial shape. 

Preliminary STAR data on azimuthal oscillations of the HBT radii [HIEI confirm 
the oscillation pattern in the left panel of Fig. [3] but disagree with the one shown in 
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Figure 3. Oscillations of the HBT radii for different transverse pair-momenta in 
RHIC collisions (left) and for the source elongated into the reaction plane (right) 
1111 . The thin circled lines in the right panel show the geometric contributions to 
the HBT radius parameters. 



the right panel. This suggests that hydrodynamics gets the overall duration of the 
expansion stage about right, since significantly longer times until freeze-out should 
cause the source to reverse the sign of its spatial deformation. More quantitative 
statements have to await a resolution of the above-mentioned "HBT puzzle" in central 
collisions. 

The significance of the i^x-dependence of the oscillation amplitudes for the HBT 
radius parameters is still largely unexplored. From the analysis of central collisions 
|2U| it is known that the ^-averaged transverse HBT radii (in particular R ) exhibit 
stronger dependence for sources with a sharp surface (such as the hydrodynamic 
one obtained from the Cooper-Frye prescription) than for Gaussian sources. One 
might expect similarly strong differences for the fCx-dependence of the oscillation am- 
plitudes in non-central collisions. 
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